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Analyzing correlation functions of charmonia at finite temperature (T) on 32 3 x (32-96) 
anisotropic lattices by the maximum entropy method (MEM), we find that J/ip and rj c survive 
as distinct resonances in the plasma even up to T ~ f.6T c and that they eventually dissociate 
between 1.6T C and 1.9T C (T c is the critical temperature of deconfinement). This suggests that the 
deconfined plasma is non-perturbative enough to hold heavy-quark bound states. The importance 
of having sufficient number of temporal data points in MEM analyses is also emphasized. 



PACS numbers: 12.38.Gc,12.38.Mh,25.75.Nq,14.40.Gx 

Whether hadrons survive even in the deconfined quark- 
gluon plasma is one of the key questions in quantum chro- 
modynamics (QCD). This problem was first examined in 
Q and in different contexts. In the former, it was 
shown that collective q-q excitations with a low mass and 
a narrow width in the tt-u channels exist even above T c 
(the critical temperature) from analyses of the spectral 
functions in the Nambu-Jona-Lasinio model. The fate of 
heavy mesons such as J/i/j in the deconfined plasma was 
also investigated in a phenomenological potential picture 
taking into account the Debye screening |3j. In general, 
there is no a priori reason to believe that the dissociation 
of bound states should take place exactly at the phase 
transition point |4j. 

Experimentally, measurements of dileptons (dipho- 
tons) in heavy ion collisions may provide a clue to the 
properties of vector (pseudo-scalar) mesons in hot/dense 
matter. Indeed, data from CERN SPS indicate anoma- 
lies in the dilepton spectra relevant to p and J/ij). Also, 
RHIC is going to produce ample data of dileptons in a 
few years j5j- 

From the theoretical point of view, the spectral func- 
tion (SPF) at finite temperature T, which has all the in- 
formation of in-medium hadron properties, is a key quan- 
tity to be studied. Recently, the present authors have 
shown [|| that the first-principle lattice QCD simulation 
of SPFs is possible by utilizing the maximum entropy 
method (MEM). We have also formulated the basic con- 
cepts and applications of MEM on the lattice at T = 
and TV in 0. The purpose of this Letter is to report 
our latest results of SPFs of low lying charmonia (J/tp 
and rjc) above T c extracted from quenched QCD simu- 
lations on anisotropic lattices. We find that J/i/j and 
i] c survive even up to T ~ 1.6T C and disappear between 
1.6T C and l.9T c . This together with our previous results 
for ss mesons on the same lattice Q suggests that the 
system is still strongly interacting above T c . 

To draw the above conclusion with a firm ground, we 
put special emphases on [I] the MEM error analysis of 
the resultant SPFs and [II] the sensitivity of the SPFs to 



Adata (the number of the temporal data points adopted in 
MEM). These tests are crucial to prevent fake generation 
and/or smearing of thepeaks and must be always carried 
out as emphasized in 0H| • 

Let us first summarize the basic formulation of MEM 
applied to lattice data at finite T |2|. We consider the 
Euclidean correlation function of the local interpolating 
operator Ji (t, x) = cijic with I — 1,2,3 for J/ip and / = 5 
for rj c . The spectral decomposition of the correlator in 
the imaginary time < r < l/T reads 

D(t) = y"(J ; (r,x)j/(0))d 3 x = ^ K(r,u)A((j)du, (1) 

where w is a real frequency and A(u>) is the spectral func- 
tion. The sum over I = 1, 2, 3 is taken for J/ip. K(t, lj) is 



the integral kernel, K(t, uj) 



-(l/T-7> 



)/(l 



). For simplicity, we take the three momentum of 
the correlation function to be zero. 

Monte Carlo simulations provide Dfa) with statistical 
error on a discrete set of temporal points Tj. Although 
there exist infinitely many A(ui) which give the same 
D(ti) through the integral Eq.QJ, MEM provides a way 
to select a unique A(ui) without introducing parametriza- 
tion of SPF and to give its statistical significance on the 
basis of the Bayes' theorem. The most probable A(ui) 
given lattice data D is obtained by maximizing the condi- 
tional probability P[A|D] cx e aS ~ L , where L is the stan- 
dard likelihood function and S is the Shannon- Jaynes 
entropy: 



A(w) - m(u}) - A(u>) lof 



(A{u 



\m(u) 



duj. (2) 



Here a is a parameter dictating the relative weight of 
S and L. The statistical significance (error) of the 
resultant A(u>) is estimated by the second variation, 
(S/SA) 2 P[A\D]. The default model m in Eq.©, the 
first estimate of which is obtained from the perturba- 
tive QCD calculation of A(u> 1 GeV), may be cho- 
sen so that the MEM errors become minimum. The fi- 
nal result A(ui) is given by the weighted average over 
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a: A(uj) = J A a (ui) P[a\Dm] da, where A a {u>) is ob- 
tained by minimizing P[A|-D] for a fixed a. The condi- 
tional probability P[a|Dm] can be calculated by using 
the Bayes' theorem and the lattice data. Thus a is even- 
tually integrated out and does not appear in the final 
result. 

MEM has been successfully applied to the lattice QCD 
data at T — to extract theparameters of the ground 
and excited state hadrons [6J, IsL llfj • On the other hand, 
applications to T ^ system have been known to be 
a big challenge [g, although there are some studies at 
high T ^lj. The difficulty originates from the fact that 
the temporal lattice size L T is restricted as L T — l/T = 
N T a T , where a T (N T ) is the temporal lattice spacing (the 
number of the temporal lattice sites). Because of this, 
it becomes more difficult to keep enough A^ata to obtain 
reliable SPFs as T increases. In other words, simulations 
up to a few times T c with A^ata as large as 30 require 
fine lattices at least in the temporal direction (anisotropic 
lattices). 

On the basis of the above observation, we have car- 
ried out quenched simulations with (3 = 7.0 on 32 3 x N T 
anisotropic lattice. The renormalized anisotropy is £ = 
a a /a T — 4.0 with a a being the spatial lattice spacing. We 
take the naive plaquette gauge action and the standard 
Wilson quark action. The corresponding bare anisot ropy 
£0 = 3.5 is determined from the data given in 0, Il3| . 
The fermion anisotropy jp = n T /n a with k„ (k t ) being 
the spatial (temporal) hopping parameter is determined 
by comparing the temporal and spatial effective masses 
of the pseudoscalar and vector mesons on a 32 2 x 48 x 128 
lattice. a T — a a /A — 9.75 x 10 -3 fm is determined from 
the p meson mass in the chiral limit |2(ij |. The physi- 
cal lattice size in the spatial direction L a = 1.25 fm can 
well accommodate J/4 or rj c , whose root mean square 
radius is about 0.5 fm. N T , the corresponding T/T c , and 
the number of gauge configurations N gSuUge are summa- 
rized in TABLE I. Calculations close to T c (T/T c = 1.04 
and 0.93) are also currently undertaken. Gauge configu- 
rations are generated by the pseudo heatbath and over- 
relaxation algorithms with a ratio 1:4. Initially, the gauge 
field is thermalized with 10000 sweeps and, then each 
configuration is separated by 1000 sweeps. 



TABLE I: Temporal lattice size, the corresponding tempera- 
ture, and the number of gauge configurations in our simula- 
tion. 





32 


40 


46 


54 


96 


T/T c 


2.33 


1.87 


1.62 


1.38 


0.78 


Agaugc 


141 


181 


182 


150 


194 



We have measured the point-point qq correlation func- 
tions in the scalar (S), pseudo-scalar (PS), vector (V), 
and axial-vector (AV) channels for the spatial hopping 



parameters K a = 0.08285, 0.0850, 0.0853, and 0.08545. In 
this Letter, we focus our attention on the heaviest quark 
mass (n a = 0.08285 with 7^? = 3.476), and study the V 
{J/4 1 ) and PS {r) c ) channels. The masses determined on 
the T = lattice (32 2 x 48 x 128) are mj/^ = 3.10 GeV 
and m Vc = 3.03 GeV. 

In all the figures below, we adopt the continuum kernel, 
K(t, u), defined after Eq.QJ. Use of the lattice kernel 
0, QjJ does not lead to any appreciable difference. This 
is because our lattice spacing is quite small, a T ~ 0.01 
fm [2l]]. Following [tJ, we define dimensionless SPFs: 
A{ui) = Lo 2 p{ijj) for T] c and A{ui) = 3uj 2 p{uj) for J/4- If 
the temporal distance between the source and the sink is 
closer than £a T , lattice artifact due to anisotropy would 
appear in the SPFs for lu > 7r/£a r . To avoid this, we 
exclude the six points near the edge (tj = 1,2,3 and 
N T — 3, N T — 2, N T — 1) and adopt the points r, = 4, 5, • • • 
and N T — 4, N T — 5, • • • until we reach the total number 
of points 7V data (< N T -7). 

The current operator on anisotropic lattice jj jAT {x) 
and that in the continuum Jp ON (a;) are related by 
jLAT ^/ aTa 5 jCONy (2Zii/Ki-K a ), where Zi is a nonper- 
turbative renormalization constant. The default model 
m{bj) extracted from the perturbative evaluation of 
A{oj 1 GeV) also depends on Zi 0. We vary Zj from 
1 (the value in the weak coupling limit) and the known 
value at j3 = 6.0 on the isotropic lattice The cen- 

tral values of m thus obtained are m(w)/(3w 2 ) = 0.40 and 
m{uj)/oj 2 = 1.15 for the V and PS channels, respectively. 
We use these values throughout the following figures. We 
have checked that our conclusions are not modified within 
the variation of Zi mentioned above. 

Shown in FigHJare p{u)s for J/4 at T/T c = 0.78, 1.38, 
and 1.62 (FigOJa)) and those at T/T c = 1.87 and 2.33 
(FiglUb)) (The corresponding Adata used in these figures 
are N data = 89, 40, 34, 33, and 25 from low T to high T). 
If the deconfined plasma were composed of almost free 
quarks and gluons, SPFs would show a smooth structure 
with no pronounced peaks above the qq threshold. To the 
contrary, we find a sharp peak near the zero temperature 
mass even up to T ~ 1.6T C as shown in FigQJa), while 
the peak disappears at T ~ 1.9T C as shown in Fig^Jb). 
The width of the first peak in Fig^a) partly reflects the 
unphysical broadening due to the statistics of the lat- 
tice data and partly reflects possible physical broadening 
at finite T. At the moment, the former width of a few 
hundred MeV seems to dominate and we are not able to 
draw definite conclusions on the thermal mass shift and 
broadening. 

The second and third peaks in Fig^a) may be re- 
lated to the fermion doublers as first pointed out for light 
mesons at T = . This must be checked by studying 
whether the peak position scales as \/a a by varying a a , 
which remains to be a future problem. Our error anal- 
ysis shows that the peak around u = in Fig^b) is 
not significant in the present statistics. The qualitative 
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change of the spectral structure between T/T c = 1.62 
and 1.87, and other features seen in the J/ip channel are 
also observed in the r] c channel as shown in Fig|2J 

Let us now evaluate the reliability of the existence (ab- 
sence) of the sharp peak at T/T c — 1.62(1.87) by the two 
tests [I] and [II] mentioned before. The first test is the 
error analysis of the peak. Shown in Fig|3|are the SPFs 
for J/ip at T/T c = 1.62 and 1.87 with MEM error bars 
(The frequency interval over which the SPF is averaged 
is characterized by the horizontal position and extension 
of the bars, while the mean value and the la uncertainty 
of the integrated strength within the interval are charac- 
terized by the heights of the bars). The sharp peak at 
T = 1.62T C is statistically significant, and the absence of 
the peak at the same position at T = 1.87T C is also sta- 
tistically significant. The same features are also observed 
for rj c . 

The second test is the A^ata dependence of the SPFs. 
Shown in FigQJa) is a comparison of the SPF obtained 
with -/Vdata = 34 and that with A/^ata — 39 for the same 
temperature T — 1.62T C (N T = 46). The two curves are 
almost identical with each other (Note that the maximum 
number of temporal data available in this case is 46 — 7 = 
39). FigEJb) shows SPFs obtained with 7V data = 26 and 
33 for a higher temperature T = 1.87T C (N T = 40). The 
maximum number of data points available is 40 — 7 = 33 
in this case. Again the two curves are almost identical. 
Therefore, the qualitative change of the SPF between 
T = 1.62T C and 1.87T C is the real thermal effect and not 



caused by the artifact of the insufficient number of data 
points. 

Here we make a brief comment on related works ai ming 
at studying the charmonia above T c using MEM llg • 
In contrast to our large temporal grids (N T = 46 for 
T = l.62T r ), N T in these papers are 3-4 times smaller 
(N T = 17 in (HI and N T = 12 in 16] at T = 1.62T C ). For 
such small N T with point-point correlation, the reliability 
tests [I] and [II] need to be done before drawing physics 
conclusions. There is indeed evidence that N T = 17 fails 
the test [II] as shown in |l5| . 

In summary, we have extracted the spectral functions 
of J/ip and rjc in the deconfined plasma using lattice 
Monte Carlo data and the maximum entropy method. 
In the quenched approximation, the number of temporal 
sites N T is taken as large as 46 and 40 for T/T c = 1.62 
and 1.87, respectively. Careful analyses of the MEM er- 
rors and the A^ata dependence of the results are carried 
out. It is found that there are distinct resonances up to 
T ~ 1.6T C and they disappear between 1.6T C and 1.9TU. 
This together with our previous results on ss mesons |8j 
on the same lattice indicates that the quark-gluon plasma 
is still strongly interacting above T c so that it can develop 
hadronic resonances. Whether these resonances found in 
quenched simulation survive in full simulation with dy- 
namical quarks is an interesting future problem. It is 
also important to unravel the nature of such resonances 
by studying their spatial structure on the lattice or 
by using non-relativistic potential approaches [lsj . 
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FIG. 1: Spectral functions for J/ip (a) for T/T c = 0.78, 1.38, 
and 1.62 (b) for T/T c = 1.87 and 2.33. 



FIG. 2: Spectral functions for r) c (a) for T/T c = 0.78, 1.38, 
and 1.62 (b) for T/T c = 1.87 and 2.33. 
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FIG. 3: Spectral functions for J/ip with MEM errors (a) for 
T/T c = 1.62 (b) T/T c = 1.87. 
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